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INTRODUCTION 


The Htiidy of inmteady turbulent boundary- layer flows has become very 
im|)ortant during the past several years. Physical experiments In this area 
are difficult to perform because of the complex meclianisms and sophisticated 
instrumentation that are re<iulred. As a result, accurate measurements are 
scarce. Most of the research el fort has been expended on developing c«imputa- 
tlunal methods for modeling the related equations, and several numerical 
experiments have been performed by various authors. However, these numerical 
studies have seldom addressed the same problem; as a result, little detailed 
comparison of results has been possible. In the rare case when a physical 
experiment has been available, the results of the various methods have 
differed widely, even on so simple a case as oscillating flow on a flat plate 
(fig. I). Tliere renuiin many unanswered questions, and a more coordinated 
effort is needed so that the effect of numerics can bJ Isolated from the 
effects of turbulence modeling, as well as from ?• xper Imcntal inaccuracy. 

Five numer'eal experiments are proposed liercln as standard cases to be 
studied by all who attempt the analysis of unsteady turbulent boundary- layer 
behavior. Iliese five problems do not encompass all of the fluid mechanics 
associated with unsteady boundary- layer characteristics; however, they are 
suggested as a common starting point for all methods, so that differences 
between the results of varioius mc'thods can be studied in a completely defined 
environment, li/lien the initial and boundary conditions are fully defined, the 
variances In the results can be studied with a much reduced range of unknown 
parameters to evaluate. 

TIk* test cases are thus presented as an Idealized set of experiments 
analogous to those of the 1968 Stanford conference (ref. 2), except that here 
the experimental data are mathematically 'eflned. llie various methods can 
then be compared on common grounds, and the results Interpreted in a known 
environment. Tlie test cases have been designed to increase progressively in 
difficulty and thus test systematically the capabilities of proposed numerical 
methods. Tliey include steady flow on a flat plate, steady adverse-pressure- 
gradient flow, oscillatory flow on a flat plate, unsteady adverse-pressure- 
giadlent flow, and ultimately, unsteady reversed flow. 


STtJVDY n.()W US INC TIMK-KKLAXATION 


CaHc 1: Zero PreHsiire-Ciratl lent 


Steady turbulent boundary- layer flow on a flat plate liaa been chosen as 
the first test case because It allows evaluation of candidate computer pro- 
grams In a well-known envl ronim-nt which has been accurately mtideled by many 
8teady-fli>w methods (ref. 2). In addition, the results of these calculations 
are to be used as starting conditions for fltjw Cases 3-5. Tlie chosen test 
case Is the Uleghardt flat-plate flow No. 1400 as evaluated In reference 2. 

In order to permit programs to develop properly the requisite Initial profiles 
for pertinent parameters, computation imiy begin ahead of the specified start- 
ing station and proceed to the starting location while permitting any tran- 
sient effects of the particular technique to dissipate before computing the 
required run. The starting condition of Rn ■ 5000 (x ■ 1.24 m) has been 
selected to reduce the possible contribution of low Reynolds number effects to 
the problem; the calculation Is to progress to x • 4.69 m (Rej. ■ 10^). 


Required itqmts- Match or Input the experimental u/U^. vs y/6* at 
■ 1.24 m; Impose Ui.(x) of reference 2 as external velocity condition. 


Rotiuired prooedure- Using the velocity profile u/U^ vs y/6* at 
X ■ 1.24 m as initial conditions at all subsequent x calculation stations, 
compute the flow as a time-dependent relaxation to steady-state conditions. 
Indicate the time required for the 6* to attain 9931 of Its computed steady- 
state value at x “ 4.69 m. Tills time should be presented In terms of 
reference lengths traveled by a particle moving at free-stream velocity 
(1. ■ 3.45 ml. Since this Is also a test of the program's capability to 
adjust to variations in initial conditions, an Indication of the size and 
number of computational steps should be presented. Convergence must be dem- 
onstrated In the variation of 6* vs t for a time equivalent to 0.25 (l./lD 
after 6* has reached 99%. 


Required cfompari eoria~ Comparisons with experimental data of reference 2 
are required for the velocity profiles at x “ 1.24 m and 4.69 m; note that 
all y profiles are to be scaled by displacement thickness. In addition, 
the X distributions of wall shear stress, displacement thickness, and shape 
factor are required. Optional: Present the variation of shear stress 

vs (y/6*) at the initial and final stations. 


Case 2: Adverse-Pressure-Gradient Flow 

Tlie analysis of spatial adverse-pressure giadlents plays an Important 
role In the study of unsteady turbulent boundary-layer behavior. Therefore, 
the second test case requires demonstration that each method can properly 
represent such developments In steady flow prior to addressing the unsteady 
problem. The Bradshaw adverse-pressure-gradient Flow C (No. 33CO of ref. 2) 
Is proposed for this purpe-sc. 
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h't’ifuii'thl htf'utif- Matrli or input tin* oxpor Imontul u/U«. vh y/6* at 
X ■ m (2.0 ft); impotif Uy(x), .ind DU^.(x) as kIvoii In rofiTonro 2 .ih 

oxtornnl volurity and dll^./dx condltlonK. 

H>'quii'i'd pri'i'i - Uko tin* voloclt*' profilt* u/U^. vh y/fi* at 
X ■ ().(>! m AH lilt* Initial rondltl(»n for all computation Htations from 
X • O.bl m to X • 2.1) m (7.0 ft). KollowinR the method of CaHe 1, compute 
the flow stiowiiiK the time required for 6* at x • 2.1 ) m to reach 99Z of 
the computed Hteady-Htate value. if a program Ih not capable of accepting 
thin HtartinK condition, the immerical experlmiuU can be ntarted usiiiK the 
tabulated experimental profile data of reference 2. In either cane, Hhow 
an a function of time for 0.25(l./l>„) past tlie time noted above to demonstrate 
converjience . 

Hequt t'i'd i^ona- (aimpar Ihimih with experimental data of referenrt* 2 

are required for the velocity profiles plotted vs y/i* at x ■ 0.61 ra and 
X » 2.31 m. The x dlHtributlons of wall she.ir stresH, displacement thick- 
ness, and shape fatt»)r between tiie beKinnlnK and ending stationa are required. 
Optional: I'resent tlie variation of the shear Htress plotted vs (y/fi*) at the 

Initial and final stations. 


unsti:ady flows 


Case ): Oscillating Flow on a Flat Plate 

Tlie oscillating flow over a stationary flat plate offers the simplest 
extension of turbulent boundary-layer methods to unsteady flow, both physically 
and m;itl)emat leal ly , because the turbulent structure and scaling laws arc not 
complicated by the additional effects of spatial pressure gradients. It is 
also a case for which some experimental data already exist (refs. J and 4). 

Kyen in tills case there is a gre.it lilverslty of results, as shown in figure 1. 
However, the differences may possibly be attributed to different interpreta- 
tions of the initial and boundarv conditions, a difficulty that will be 
alleviated by the use of the projiosed standard-case conditions. 

Since the steady flow on a flat plate Is a universal starting place for 
turbulent boundary-layer theories, the well-documented Wieghardt flow of 
Case 1 has been selected as the initial conditions for the proposed oscil- 
lating flow test. Again, the starting condition Is chosen at an R() such 
that no low Reynolds number effects will be present. Note: the reference 

length for reduced fr'»quency computation corresponds to the length of the 
unsteady portion of tne flow rather than the dist.ince from x ■ 0. This 
length has been chosen because the flow is considered to be quast-steady up to 
X “ 1.24 m; therefore, this initial distance will not contribute to unsteady 
effects. 
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Hcquli\'d inputiJ- Tlu' ri*«iul r«-d InputH ;irf: 

1. t 5 0, X - 1.24 m: WloKhardt H.«t-platf prof I It- as dflt-rml lud for 

Cast* 1, with R(i • 4860, modified to forroHpond to hual idm-ve lot 1 1 Ion such 
that R() - 4Hh0 (U^./L’«) ° 

2. X i 1.24, t • 0: roHultK of computation of Wlcuhardt flat-plate 

flow as determined In Case 1 for x • 1.24 m t«> x ■ 4.69 m 

1. Ed^e conditions: Uj^,/lU ■ tl.O + A sin wt) 

4. Test conditions: u ■ O.S, 1.0, 5.0, 10.0, where u) ■ litl./L'u,, 

U«, • 33 m/s, 1. • 3.45 m; A ■ 0.1 for all f retinenc les (Optional: A ■ 0.5) 

Requiivd Using the specified Inputs, compute three cycles of 

oscillation ft>r each fretiuency. All output data are to he based on the third 
cycle of computed results, with Infornvitlon available at 30 incremt*nts 
through this cycle. 

Hi^quityd output-o- Tlte required o«itputs are: 

1. Fourier analysis of the velocity profile at x ■ 4.29 m, presenting 
the mean and first harmonic as a function of y/6* for each frequency (see 
fig. 2). When available, present the shear stress profile at x - 4.29 m In 
the same m.inner for each frequency. 

2. Plots of the wall shear stress, displacement thickness, and shape 
factor as functions of uit for the Last cycle, showing all frequencies 
computed . 


3. Plots and tabulation of phase and amplitude of wall shear stress, 
displacement thickness, and shape factor as functions of w(x - 1.24)U.„ up 
to ujl./Uuo for each frequency. Optional: Plot of the phase data as a func- 

tion of (x - 1.24) showing the development of the ph.ase angle along ^1 ( late 
up to < 1 ) 1 , /lU. 


Case 4: Unsteady Adverse-Pressure-tJradient (Unsteady Howarth Flow) 

Although some experiments are in progress (cf. ref. 5), no fully docu- 
mented experimental data presently exist tliat demonstrate the effects of 
strong, unsteady .-idverse-pressure-grad lent s on turbulent boundary layers. 
Previous analyses have considered flows with complicated combinations of 
spatial and tempor.al pressure gradients (refs. 6 and 7), making It difficult 
to assess the effect of unsteadiness on specific boundary- layer properties. 

In contrast, the unsteady Howarth flow chosen here starts from a well- 
established steady flat-plate flow, on which a linear deceleration of U^. is 
Imposed at t “ 0. After t ■ 0, the value of dU^/dx Is Increased linearly 
with tlmi-. This case has been selected In order to simplify the relation 
between spatial and temporal pressure gradients while demonstrating the 
Important features of unsteady turbulent boundary layers experiencing adverse 
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preMHuri! gradients. The test cnni will only be culculaled to the point of 
flow reversal at the Kurtace, so th.it no difficulty is introduced by ill- 
defined dovmst ream boundary conditions. 

A primary effect of unsteadiness on boundary- layer flows experiencing 
adverse pressure-gradients is to delay and even to reduce the importance of 
the strong singularity normally occurring at flow reversal (ref. H) in steady 
flows. In the unsteady environment, this singular behavior does not always 
appear and the flow near the wall can stagnate witliout strong changes in the 
character of the integral parameters. (This is illustrated by the solid 
curves in figure i for the conditions described below.) 

Hcquii’Pti inputit- The required inputs are: 

1. t i 0, X ■ 1.2A m: Wleghardt flat-plate profile ns determined for 

Case 1 at x ■ 1.24 m (R^^ * 4860), maintained as a quasi-steady input profile 
for all subsequent time 

2. X > 1.24, t “ 0: results of computation of Wieghardt flat-plate 
flow as determined in Case 1 from x ■> 1.24 m to x ■ 4.6^ m 

3. Edge conditions: Ug/U.B ■ (1.0 - AXt), 1.24 i x i 4.69 m, 

X - (x - 1.24)/L 

4. Test conditions: Ua, ■ 33 m/s; L ■■ 3.45 m; A ■ 2.4 sec”' 

It may be noted that for the values of the constants chosen, a fluid 
element in the free stream will have moved approximately 4L by the time the 
local edge-velocity at the final x station reaches stagnation. 

Requin^d pvooeduiSie- The flow is to be assumed quasi-steady up to 
x “ 1.24 m; the velocity distribution is then imposed as a linear function of 
X and t until the wall-shear stress equals zero at the final station. Two 
conditions are to be studied: (1) compute the flow until the wall shear 

stress reaches zero at the final station under the Influence of the fully 
unste.ady free-stream flow; and (2) calculate the unsteady turbulent boundary 
layer that develops when the l)^.(x) distribution that results at the com- 
pletion of condition (1) is suddenly imposed on the flat plate flow of Case 1 
Compute the flow until the X location of the flow reversal point stabilizes 
Indicate the time at which separation is considered to have occurred. 

Requirc'd outpuLa- The required outputs are: 

1. Velocity profiles vs y/6* at x ■ 2.29, 2.89, 3.79, and 4.69 m 

(X ■ 0.3, 0.48, 0.74, and l.O) at the completion of condition (1), and the 
final result for applicable X values for condition (2). Optional: Present 

the shear stress vs y/6* for each of the above conditions. 

2. Present the variation of wall shear stress, displacement thlc'tness, 
and shape factor as a function of x from 1.24 to 4.69 m for the two con- 
ditions cited (sec fig. 3). 
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C.iHi* 5: UnHtfady Flow-Ki*viT8al Without Separation (Constant- 

Foeus Hyperbola Unateady Turbulent Bubble) 

A iw>dlfled form of the unsteady turbulent bubble originally proposed in 
reference 9 has been selected as the test case for evaluating tiie use of 
unsteady boundary- la/er methtxis for analysis of flows where sign! f leant flow 
reversal appears before separation occurs. Tlie bubble overcomes one of the 
major difficulties in studying reversed flows in steady boundary layers; 
namely, that tijere is seldom a realistic way to model the incoming flow from 
the down-stream boundary of the separated region. Tlie constant-focus hyper- 
bola has been selected becr.use the discontinuity at the velocity minimum of 
tlie original case of reference 9 has been found to contribute significantly to 
flow behavior in other regions of the flow field. 

Tlie calculations are to be started using a constant zero-pressure-gradlent 
fliiw (Case 1) and then the adverse-pressure-gradient is to be Introduced an a 
function of x and t. Wlien the bubble forms, the fluid elements contained 
within the unsteady turbulent bubble have all been introduced into the flow 
from the upstream boundary at an earlier time. Therefore, the history of all 
the fluid elements Is known throughout the calculation, and the fluid elements 
can be traced as they are brought to rest and reverse within the bubble. The 
fluid elements vihin the reversed-f low region retain their gross turbulence 
characteristics, and therefore an eddy-viscosity as well as more sophisticat''d 
modeling may apply. 

Tlie selected problem is one which encompasses a significant region of 
reversed flow where integral parameters have not been grossly distorted, as 
well as a region of reversed flow where sucli distortion does occur. The shape 
factor is of particular interest, because it seems to indicate a demarcation 
between regions of simple reversed flow and reversed flow where the momentun 
thickness has significantly Increased, thus causing a rapid drop in shape 
factor (fig. 4). Tills test case is suggested as a case that may help deter- 
mine the limitations of boundary- layer theory in the study of unsteady 
turbulent boundary layers. 

Required inputs.- Tlie required inputs are: 

1. t > 0, X ■ 1.24 m; Wieghardt flat-plate profile as determined for 
Case 1 at x ■ 1.24 m (Rq ■ 48b0) , maintained as a quasi-steady input profile 
for al 1 subsequent time 

2. X > 1.24, t ■ 0: results of computation of Wieghardt flat-plate 
flow as determined in Case 1 from x ■ 1.24 m to x ■ 4.69 m 

3. Edge conditions: Ug/Uo, ■ 1.0, x < 1.24 m, 

Ue/U„ - 1.0 + I a 2 + (Bt)^(X - Xo)^! 1/2 - ( a 2 + (BXot)2) 1/2, X - (x - 1.24)/L, 
1.24 m s X i 4.69 m 

4. Test conditions: Ua. ® 33 m/s; L “ 3.45 m; A ■ 0.05; B ■ 3.4 sec“l; 

Xo - 0.7 
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hi'<iuivod pmoeduivc AMSiinH* flow 1m quaHl-ntundy up to x ■ 1.24 m, then 
impoHc the c««nBtant- focim hyperbola velocity dlBtrlbutlon as a function of 
X .ind t. The final time will be determined by the performance of each 
individual proKram. If reaulta indicate that the developmenta obnerved In 
flKure 4 occur, then compute until the plateau in ahape factor la well defined. 
Thia apec 1 f Icat ion la Intentionally vaKue, aliice thla caae la offered aa a 
teat of the lira! tat Iona of boundry- layer theory aa well aa a atnndardixed teat 
for turbulent boundary- layer ca 1 culat Iona , and no arbitrary atoppinn point can 
be accurately establlahed. 

/icv/wt JH'ci outputs The required mitputa are; 

1. Velocity profiles va y/fi* at x • 2.29, 2.89, 3.79, and 

4.ft9 m (X • 0. iO, 0.48, 0.74, and 1.0) at time intervala aaaoclated with 
0.4 1./U..,, and at last time atep calculated. Optional: I’reaent the ahear 

streaa va y/fi* for each location. 

2. Present the variations of wall ahear stress, displacement thicknesa, 
and ahape factor as functions of x at time intervals aaaoclated with 

0.2 L/U<u. Hie carpet-plot technique of figure 4 nuiy be used, or a conven- 
tional plot of the data in fixed coordinate form is acceptable. Indication 
should be tmule of the first time at which flow reversal is detected. 


CONCLUSIONS 


Five special test cases have been presented in an effort to bring all 
the various unsteady turbulent boundary- layer calculation mc-thods together, 
at least for a few common points where direct comparisons can be made. Some 
correlation with experiment in steady flow is Included to verify the basic 
performance of the various methods, but the emphasis is on analysis of the 
Importance of unsteadiness In zero pressure-gradient, adverse-presaure- 
gradient, and reversed flows. These tests are not expected to explain all 
aspects of the fluid mechanics of unsteady turbulent boundary- layer flow 
behavior; rather, they are offered as a common starting point for further 
studies. Tt is the author's hope that the use of these cases will mean that 
future numerical experiments can be better coordinated with existing know- 
ledge, that the results of these numerical experiments can be better under- 
stood, and that new methods can then be Judged based on a common standard. 
The author would appreciate receiving tabulated samples of any computations 
of these test cases; a compendium of results will be prepared for later 
publ icat ion. 
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OSCILLATING FLAT PLATE 



y 


free stream (from ref. 1). 



DETAILS OF VELOCITY AND VISCOUS PARAMETERS 



U/Ua (1 + A sin tjt) r/r 

“ u w 

6* = + Ab* sin (cut + 

H = Hq + AH^ sin (cut + 

Figure 2.- Typical flow parameters for a flat plate with an oscillating 

free stream. 


10 



WALL SHEAR STRESS DISPLACEMENT THICKNESS 


I 



1 1 


$ 

o 


« 

o 

UJ 

I- 

< 

QC 

< 

CL 

UJ 

C/5 

% 



UJ UJ 
> c/5 
UJ Z 
CC D 

□ 




12 


Figure 4.- Shap- factor as a function of time and X 

hyperbola unsteady turbuleu bubble. 
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Five numerical experiments are proposed as standard cases to be compated ' 
by all who plan to analyse unsteady turbulent boundary-layer behavior. In 
this way, differences between the results obtained by various methods can be 
compared in a completely defined environment. Tlie test cases ranRe in 
difficulty from t ime-relaxat ion study of the steady flow on a flat plate to 
the analysis of unsteady reversed flow. Initial and boundary conditions are 
fully defined for each case and representative outputs are presented. It is 
recommended that tabulated samples of computations of these test cases be 
published in a compr iditim of results. 
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